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Abstract 

We consider the asymptotic expansion of the heat kernel of a generalized Laplacian 
for t — > + and characterize the coefficients o^, k > 0, of this expansion by a natural 
intertwining property. In particular we will give a closed formula for the infinite order 
jet of these coefficients on the diagonal in terms of the local expressions of the powers 
of the given generalized Laplacian in normal coordinates. 
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1 Introduction 

A generalized Laplacian A acting on sections of a vector bundle £ over a Riemannian ma- 
nifold M is a second order elliptic differential operator with scalar symbol ( PGV Q. Its 



spectrum is interesting for various geometrical questions and the spectra of the more promi- 
nent generalized Laplacians are of fundamental importance in global analysis. The short 
distance behaviour of the spectrum of a generalized Laplacian can be studied using analytic 
properties of its heat operator A^ A , which solves the heat equation 



d_ 

Of 



associated to A for all sections ip of £ . The general existence result of Minakshisundaram- 
Pleijel ([p|,[ |BGY| ) shows that Kf" is an integral operator for all t > with integral kernel 



k^(x,y), which is a smooth section of the bundle with fiber Hom(£ x ,£ y ) over the point 
(x, y) of the product M x M. Except for a very limited set of examples the heat kernel is 
not known explicitly, perhaps the most important exception is the heat kernel of the flat 
Laplacian A acting on sections of a trivial vector bundle V x E over a euclidian vector space 
V of dimension n: 



,12 



1 \x-yy 
k t (x,y) = — n e 4t id E 

As can be seen from this example the heat kernel kf^(x, y) peaks on the diagonal {x = y} of 
M x M as t — > + and it turns out that the rate of its divergence on the diagonal carries a 
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tremendous amount of information about the geometry of the underlying manifold M. More 
precisely the heat kernel has an asymptotic expansion along the diagonal as t — > + of the 
form 

I dist 2 Qr,y) 

k^{x,y) -y=ne it jVa^y) (1) 

k>0 



where dist 2 (x,?/) is the square of the geodesic distance from x to y and the ak(x,y) are 
again sections of the bundle B.om(£ x ,£ y ). In a sense these coefficients cik, k > 0, describe 
how much the heat kernel of the operator A differs in divergence from the flat Laplacian A 
in various directions and this difference does not only depend on the coefficients of A, but 
involves the geometry of M in a crucial way. 

As a matter of convenience we will not use the coefficients ak{x, y) directly but will mul- 
tiply them afc(x, y) := j(x, y)a k (x, y) by the jacobian determinant j(x, y) of the differential 
of the exponential map in y at x, which is a well defined smooth function in a neighbor- 
hood of the diagonal. Off this neighborhood we may extend j(x, y) arbitrarily to a smooth 
function on M x M, because only the germs of the coefficients ak(x,y) along the diagonal 
are characterized by the asymptotic expansion ([]]) of kf* above. The main theorem of this 
article exhibits a remarkable intertwining property of these modified coefficients a^, k > 0: 

Theorem |2.1| (The Intertwining Property) 

Let A be a generalized Laplacian acting on sections ip of a vector bundle £ over M. The 
coefficients : £ x — > £ y of the asymptotic expansion of the heat kernel of A satisfy 



A^]{y) 



£ 

u=0 



(-I)' 



[A u (a^ I/ (exp y -,y)tp(exp y -)] (0) 



where a^-^exp^-, y) ip(exp y -) is considered as a section of the trivial vector bunde T y M x £ y 
over T y M and A is the flat Laplacian acting on sections of this trivial bundle. 

The rest of the article is devoted to corollaries of Theorem |2.1| . In particular we prove the 
striking fact that the intertwining property alone characterizes the infinite order jet of all 
coefficients ak(-,y) in y and on the way we will give a different proof of Polterovich's inversion 
formula ( ||P 1|| ) for the endomorphisms ak(y,y), k > 0, of the fiber £ y in a somewhat stronger 
formulation: 



Theorem [2.4| (Polterovich's Inversion Formula) 

Let A be a generalized Laplacian acting on sections of a vector bundle £ over M and let 
dist 2 (-, y) be the square of the geodesic distance to y 6 M. For all r > k > and all sections 
ip of £ the endomorphism ak(y,y) = ak(y,y) of the fiber £ y is given by the formula : 



>*V>] (y) 



£ h)' (; + J) 



:-i) 



k+l 



(k + l)\ 



A 



k+l 



dist 



21 1 
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The final result is Theorem |27?], which provides an explicit formula for the infinite order 
jets of the coefficients ak(-,y) in y in terms of the local expressions of the powers of the 
operator A in normal coordinates on M and some smooth trivialization of £ . 
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The present version of this article was written down while I was staying at the University 
of California in Berkeley, whose hospitality is gratefully acknowledged. It is only the first 
part of a projected article of three and hence will not appear separately in printed form. 

2 The Intertwining Property 

It is well known that the coefficients a^, k > 0, of the asymptotic expansion of the heat kernel 
of A are local expressions in the curvature of M and the coefficients of A (JGfl). Consequently 
we will assume throughout this section that A is defined on sections of a trivializable vector 
bundle £ over a starshaped neighborhood of the origin in some euclidean vector space V of 
dimension n. Moreover we will assume that the identity map V — > V is the exponential 
map in the origin with respect to the Riemannian metric defined by the symbol of A, a 
fortiori the symbol metric must agree with the scalar product of V in 0. The general case 
is easily reduced to this local model by taking V = T y M to be the tangent space of the 
underlying manifold at the point y under consideration and choosing the neighborhood of 
the origin and the trivializable vector bundle £ accordingly. 

The preceeding assumptions are crucial to the entire argument but single out the origin as 
a special point. Consequently all our statements below refer to the fixed target point y = 0. 
In particular the coefficients dk(x) := a,k(%, 0) will be homomorphisms a,k{x) : £ x — ► E 
from the fiber of £ over x to the fiber E of £ over the origin. Similarly we will write 
j(x) := j(x, 0) for the Jacobian determinant of the exponential map in the origin at the point 
x G V, which is the identity map V — > V by assumption. In other words the Jacobian 
determinant is characterized by dvo\ g = j(x)dx and is the square root j(x) := a/ det g{x) 
of the determinant of the symbol metric with respect to the flat scalar product. Now by 
the remarks following ( PCV| , Theorem 2.30]) we have an asymptotic expansion for the heat 
operator 



which complements the asymptotic expansion (|IJ) of its integral kernel kf~ . Using both 
asymptotic expansions we find for any smooth section ip of £ with compact support in a 
sufficiently small neighborhood of the origin 




k t (x)ifj(x) dvolg 



E 



(-*)" 



[A"V](0) 




fi>0 

e 4t (j(x)a fJi (x)ifj(x))dx 




I 
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where a(x) := j{x)dj t {x) and A is the flat Laplacian acting on sections of the trivial E— 
bundle V x E over V. Note that the asymptotic expansion of the heat kernel kf" is locally 
uniform in x according to ( |[BGV , Theorem 2.30]) and hence we may integrate it over the 
compact support of ifi to obtain an asymptotic expansion of the integral. Existence of 
asymptotic expansions implies uniqueness and so we need only sort out the different powers 
of t to prove the intertwining property of the heat kernel coefficients in the form: 

Theorem 2.1 (Intertwining Property of the Heat Kernel Coefficients) 

Let A be a generalized Laplacian acting on sections of a trivializable vector bundle £ over V . 
The coefficients a&, k > 0, of the asymptotic expansion of the heat kernel of A are sections 
of the bundle Hom(£, V x E) and so a^ij) is a well defined section ofVxE for any section 
x/j of £. The coefficients intertwine the powers of A and A: 



(-1)' 



;a^] (o) 



E 

u=0 



[A" (0^.^)1(0) 



(2) 



A direct consequence of Theorem |2.1| is that a generalized Laplacian A acting on sections 
of a vector bundle £ remembers to vanish on sections ip of £, which are harmonic in some 
smooth trivialization of £ , in the following sense: 

Remark 2.2 (Excess Vanishing for Harmonic Sections) 

Consider a section ip of£, which is a homogeneous, harmonic polynomial Qip £ Sym d V*®E 
of degree d with A(Qif)) = in some smooth trivialization $ of £. If A is a generalized 
Laplacian acting on sections of £, then [A^ifj] (0) = vanishes for all /i < d. 

Philosophically the intertwining property (fj) provides a precise geometric interpretation 
for the coefficients of the heat kernel expansion in terms of A and the flat model operator 
A for generalized Laplacians. Similar considerations should apply to other model operators 
arising e. g. in Heisenberg calculus or in other parabolic calculi. A convenient reformulation 
of this property (fj) can be given with the help of the formal power series e~ zA and e~ zA of 
differential operators on £ and V x E respectively and the generating series 



a(z 



E 

M>0 



for the heat kernel coefficients a^, k > 0. Namely the intertwining property (E3) is just another 



way to write down the equality [e z i/j](0) 



(a(z)ijj)](0) of formal power series in z 



for all sections ip, which in turn is equivalent to the commutativity of the pentagram 



Jet£° £ 



JetS°£|U 



a[z) 



Jet^E[z' 



-zA 




JeC E\ 



Elz 



cv 



4 



where ev is the evalutation at the origin. Only the infinite order jet of the sections a^, k > 0, 
of r(V, Hom(^, V X E)) in the origin can ever be sensed by jet theory of course, and so the 
commutativity of this square does only depend on the infinite order jet of the generating 
series a(z). Strikingly however the commutativity of the pentagon characterizes the infinite 
order jet of a(z) completely In fact the unique jet solution of the intertwining property 
thought of as a set of equations for the unknowns %, k > 0, is given in Theorem [2.7| below. 



At this point we digress a little bit on a very interesting property of the flat Laplacian 
A acting on functions C°°V on a euclidian vector space V. More general we can consider 
the flat Laplacian A acting on sections C°°(V, E) of the trivial i?-bundle V x E over V, but 
the auxiliary vector space E never enters into the argument directly and so we will stick to 
the case E = R. Along with the flat Laplacian A comes the operator \x\ 2 of multiplication 
by the square of the distance to the origin, in orthogonal coordinates {x^} on V these two 
operators can be written: 

By restriction they act on the subspace Sym V* C C°°V of polynomials on V and their 
commutator [A, |x| 2 ] = (—4) (N + ^) as operators on Sym V* involves the number operator 
iV : Sym V* — > Sym V*, which multiplies Sym r V* by r, shifted by half the dimension n 
of V. In other words X := \\x\ 2 and Y := |A close with H := iV + | to an s^-algebra 
of operators on Sym V*. Iterated commutators of X and Y in s^-representations can in 
general be written down using factorial polynomials [x] r := x(x — 1) . . . (x — r + 1) or 
binomial coefficients m := ^[x] r for integral r > 0. In particular the standard relation 
Y r X r if) = r! [— X] r ip for all ip with Yi/j = and Hi/j = Xip implies for the constant 
function ip := 1 G Sym V* the classical formula: 

[A r \x\ 2r ](0) = 4 r r![-f] r . (3) 

Slightly more useful for our purposes is the following derived identity: 

Lemma 2.3 For all smooth functions ij) G C°°V and all k, I > 0: 

f-1) 



Ljd A k+l ( — \x\ 2k ib' 



(0) = (-4)' 



k 



/! 



(0). 



Proof: Only a finite number of partial derivates of ijj in the origin are actually involved 
in this identity and hence we may assume that ip is a polynomial without loss of generality. 
Moreover only the homogeneous component of ip of degree 21 contributes to left and right 
hand side, which are both evidently SO ^-invariant linear functionals in ip G Sym 2l V*. 
However there is but one SO ^-invariant linear functional on Sym 2 ' V* up to scale, so that 
it is sufficient to check the identity in question, which can be rewritten as 

[A^(N 2 ^)](0) = 4 ^^ ^p + ' [A^](0), 
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for the single polynomial ip := \x\ 21 , for which it is true by the classical formula (|3|). □ 

Returning to the general case of a Laplacian A acting on sections of a vector bundle 
£ over V we recall that the generating series a(z) for the coefficients in the asymptotic 
expansion of the heat kernel intertwines the formal power series e~ zA and e~ zA of differential 
operators. Using this intertwining property (||]) together with Lemma |2l| we calculate: 



k+l 



(0 = £ 



[ -^f- A" { l -\x\ 21 a k+l ^\ 
/i! II 



(0) 



E 



(0) 



k 
At=0 



A"(a fc _ M V) 



(0). 



Combining this equation with the binomial inversion formula 

V + |^ (-\ - /S\ (r - ^ 

1=0 



2 

r — I 



I 



valid as soon as r > > we eventually arrive at the following inversion formula for the 
heat kernel coefficients a&, k > 0, of a generalized Laplacian A: 

Theorem 2.4 (Polterovich's Inversion Formula ( [pi|| , [P2|| )) 

Let A be a generalized Laplacian acting on sections of a vector bundle £ over V and let \x\ 
be the radial distance from the origin. For any section ip G T{V,£) and all r > k > we 
can compute the action of the endomorphism afc(0) of the fiber E of £ over the origin on if) 
by means of an explicit inversion formula: 



>fcV](o) 



E H)' 



r + - 
r — I 



(-l) k+l 

. (k + iy. 



A 



k+l 



l\ 



\x\ 2i ip; 



(0) 



Our normalization of the coefficients a&, k > 0, drops the factor (47r)~5 arising from 
the value of the euclidian heat kernel for the flat Laplacian A on V at the origin to have 
the intertwining property (|2|) in as simple a form as possible. In the original formulation of 
Polterovich this factor is part of a k and of course in every conceivable application this factor 
has to be reinserted by hand. 

Example 2.5 (The Local Index Theorem for the Gaufi— Bonnet Operator) 

Choose a smooth Riemannian metric g on V satisfying our standing assumptions and a 
vector bundle £ over V with a connection V. For a generalized Laplacian on £ of the form 
Aft = V*V + HF for some smooth endomorphism F of £ we evidently have 



(-l) k+l 
(k + l)\ 



V*V + hF) k+l 



-K) k _ pk {-I) 



k\ 



-A + terms of lower order in h~ 

+ differential operators of lower order 
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for all k, I > 0. Plugging this into Polterovich's Inversion formula we verify immediately 

(-h) k 

flfc(0) = — — — F k + terms of lower order in h~ 

for all k > ; which may have been anticipated from Theorem \2. J| . Remarkably this simple 
calculation together with a few general considerations concerning the order of the coefficients 
ak is already sufficient to prove the local index theorem for the Gaufl-Bonnet operator! 

In conclusion the intertwining property (Q) alone is sufficient to determine the value of 
the coefficients a^, k > 0, of the heat kernel expansion at 0. Before we proceed to show that 
in fact the jets of infinite order jet[f a/-, k > 0, are determined by (0) we want to make a 
few general remarks concerning our guiding philosphy in the calculations to come. We want 
to avoid formulas involving compositions of differential operators, because such formulas 
are hardly if ever of any use in explicit calculations. In favourable situations it may still 
be possible to calculate the values evD r , r > 0, of powers of a differential operator T> in 
the origin without knowing the partial derivatives of their coefficients. Note that the usual 
arguments of symbol calculus become meaningless if we have no control over the partial 
derivatives of the differential operators in a singular point like the origin. 

In general the value of a differential operator T) acting on sections of £ at the origin 
will be an element of Hom( Jet|f S, E). For the trivial E'-bundle V x E however we may 
identify JetjftV x E) with the formal power series completion of Sym V* ® E and hence 
Hom( Jet^°(V x E),E) with Sym V ® Endi? in the usual way. The scalar product of the 
euclidian vector space V extends to a scalar product on Sym V* defined via Gram's perma- 
nent and characterized by (e a , e /3 )s ym y* = e^'^ for all a, (3 E V* with a slight abuse of 
notation. Alternatively we may choose orthonormal coordinates {x^} on V and write down 
the scalar product directly 



r. 

r>0 



gr 

-if; 



dXf J-1 . . . dx^ r 
((ev <g>ev) o e (v,v> ) O ® 



(0) 



dXf J-1 . . . dx^ r 



(0) 



where (V, V) is the bidifferential operator if> <g> ip \ — > -^-ip ® ^f - ^ an< ^ ev ^ s ^ ne 
evaluation at as before. Written in this form it is clear that the scalar product extends to 
the formal power series completion of Sym V* or even to smooth functions provided that the 
defining sum converges. Moreover the musical isomorphism ft : Sym V — > Sym V* of this 
scalar product is the natural extension of the musical isomorphism of V and extends to the 
formal power series completion of Sym V, too. The value of a differential operator T> acting 
on sections of V x E has an image (ev T>y under this musical isomorphism, which is the 
unique polynomial on V satisfying [Dip] (0) = ( (ev £>*), ip )s ym v* for all smooth sections ip 
of V x E. 

The scalar product on Sym V* will enter the formulas through the operator (V, V) intro- 
duced above, which in turn makes its appearance via Green's identity for the flat Laplacian 
A and the multiplication map m 

(Aom)('0®'0) = (m o (A® id - 2(V, V) + id® A)) (ip®4>) 
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which features three commuting operators A <g> id, id <8> A and (V, V }. Hence we are free 
to put these operators in arbitrary order upon exponentiation. For the ensuing calculations 
we need to choose a trivialization $ of £, which we think of as a family of homomorphisms 
: £ x — > E with $(0) = id^, so that Qip is a smooth section of V x E for every section 
i\) of £. With this in mind we find 

[e- zA (a(z)i))] (0) = (ev o e~ zA o m) (a(z)®- 1 ® <S>ip) 

= ((ev®ev) oe 2z(v ' v) ) ( e^K^ 1 ) <g> e~ zA (^) ) 

= ( (2zf (e' zA (a(z)®- 1 ), (e~ zA (<^)) ) Sym y* 

= < e z ^ 2 (2z) N (e~ zA (a(^)$- 1 ), ^ ) Sym y* 

while [e~ 2A -0] (0) = ((eve~ 2<I>A<I> $ - 0)s y mV* by definition. As the scalar product (, ) 
is non-degenerate on Sym V* and ip can be chosen arbitrarily we conclude: 

(eve- zMrl )» = e z ^ 2 (2zf (e~ zA (a(z)®- 1 )) . (4) 

Evidently the operators e z ^ 2 and e~ zA are invertible with inverses e -2 ' 1 ' 2 and e zA respec- 
tively, thus we can solve equation (f|) uniquely for the infinite order jet of a(z) leading to: 

jetg (a(^)$- 1 ) = e zA (2z)- N e^ 2 (ev e-^' 1 f 

= e zA (2zr N (eve~ z " A ^ 1 e zA f. (5) 

I 1 2 

Concerning the last step we remark that the polynomial e~ z < x < corresponds to the differential 
operator e zA with constant coefficients. Hence it is possible to replace e~ z ^ (ev e _2 * A<1, )tt 
by (ev e~ z ® A ® e zA )' ) , but not by (ev e 2A e - 2 * A,I> )tt ; because the latter would involve par- 
tial derivatives of the coefficients of e _z * A * in the origin. Nevertheless we cheated a little 
bit, because while certainly injective the operator (2z) is definitely not surjective. It is 
quite remarkable in itself that the value of the differential operator e - z ^^ e zA i n the ori- 
gin lies in the image of (2z) N by equation (|]), because its coefficients for the different powers 
z r , r > 0, of z must have order less than or equal to r to have this true! Hence the inverse 
(2z)~ N is well defined and refers to the unique preimage of (ev e _2<I>A * e zA )" under [2z) N . 

Definition 2.6 (The Difference Operator of a Generalized Laplacian) 

Consider a generalized Laplacian A and the flat Laplacian A acting on sections of £ and 
the trivial vector bundle V x E over V respectively. Associated to A and some smooth 
trivialization $ of £ is the formal power series 

R ( z ) := e - 2 (*A*-) e *A = 1 + zRi + z 2 R2 + z 3 Jh + ... 

of differential operators on V x E called the difference operator. The values ev JX r of its 
coefficients have order less than or equal to r, so that ev H(z) lies in the image of (2z) N . 

At least in principle the recursion relation -^H(z) = H(z)A — ($A$ _1 )II(,z) or equiva- 
lently rJX r = H^xA — ($A$ _1 )H r _ 1 for r > 1 together with H = 1 allows us to calculate 
JX(z) to arbitrary order in z. However the composition ($A$ _1 )H r _i requires knowledge 
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of the full operator H r _i and not only of its value at the origin. Note that the recursion 
relation in itself does not imply that the values of the H r , r > 0, have order less than or 
equal to r, because the origin will be a singular point for the full operator H(z) in general. 
Reformulating equation (|5]) according to the definition of H(z) we arrive at the following 
theorem, which in a sense is the converse of Theorem |2.1| : 

Theorem 2.7 (Characterization of the Heat Kernel Coefficients) 

The infinite order jet of the generating series a(z) for the heat kernel coefficients of a gener- 
alized Laplacian A acting on sections of a vector bundle £ over V is uniquely characterized 
by its intertwining property [e~ zA vp] (0) = \e~ zA (a(z) iJj)~\ (0). In fact the infinite order jet 
of any solution a{z) to this equation is given in terms of some smooth trivialization $ of £ 
and the corresponding difference operator H(z) by: 

jet~(a(^)$- 1 ) = e zA (2z)~ N (ev R(z))K 



The essence of Theorem [2.7| is not particularly evident in the abstract formulation above 
and so we want to close this section expanding this formula for the jets of the heat kernel 
coefficients into powers of z. This exercise will show in particular that it is just as difficult to 
calculate the infinite order jets of the heat kernel coefficients a^, k > 0, as it is to calculate 
the value of the powers A r , r > 0, of A in using normal coordinates and the chosen 
trivialization $ of £. In other words calculating the powers A r explicitly is the problem of 
calculating the coefficients of the heat kernel expansion en nuce! Let us start by expanding 
the value ev H(z) into powers of z and homogeneous pieces 

(evH(z)f = J2 zrR r,s, 

r>s>0 

where each JX r s e Sym s V* Cg> End E is homogeneous polynomial of order s, and employ 



Theorem 12.71 to conclude 



jetsv*) = E ^ At (^4 )S ) = E^E^ A ^ + M +2t (6) 

r>s>2t>0 ' k,l>0 t=0 

by reordering the sum and setting k := r — s+t and I :— s—2t. In particular the infinite order 
jet of the first coefficient a picks up the principal symbols evil;/, / > 0, of the operators 
^ > 0, in the origin. Now by Definition |2.6| of the difference operator H(z) we find for 
the homogeneous pieces of its value 



rI = Et-ir^M^T^tev^A^r-; 



s— 2/x 



where (ev ($A$~ 1 ) r ~ M )ij_2/i £ Sym s_2At V* ® End E is the homogeneous piece of order s — 2/z 
of the polynomial corresponding to ev ($A$~ 1 ) r-/ \ Inserting this expression into equation 
(P) above we get a closed formula for the homogeneous pieces of the infinite order jets of 
the heat kernel coefficients a&, k > 0. Though explicit this formula looks rather akward and 
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so we refrain from writing it down. In fact the main point is to understand the fantastic 
cancellation properties of the Aj,, r > 0, which make the ev H r , r > 0, operators of degree 
less than or equal to r in the first place! Without proper understanding of these cancellations 
we won't be able to take advantage of either formula to calculate the heat kernel coefficients 
explicitly. 
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